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(sj ! I Introduction 



The success of the Kobayashi-Maskawa (KM) model of CP violation [TJ in predicting 
correctly CP asymmetries in B meson decays has been recently recognized by the 
Nobel committee [2] . The large asymmetry measured in B° — > J j ipKs(L) (21 S] given 
by sin 2(3, where (3 = 4>\ = arg(— V* b V c d/V* b Vtd), proved unambiguously that the KM 
phase is the dominant source of CP violation in B decays. 

This test involves an interference between B°-B° mixing and a B° decay ampli- 
tude [5], consisting of a dominant color-suppressed b — > ccs tree amplitude and a 
small contribution from a b — > suu penguin amplitude [51 El IB]- (We use the unitarity 
of the Cabibbo-Kobayashi-Maskawa (CKM) matrix, V t * b V ts = -V* b V cs - V* b V us , 
and will not discuss order 10~ 3 effects due to CP violation in B° — B mixing [ID] and 



mixing [11].) The magnitude £ of the ratio of these two amplitudes, which 
determines the theoretical precision of this test, involves three suppression factors: 
A ratio of CKM matrix elements, \V* b V us \/\ V* b V cs \ ~ 0.02 [TJ, small Wilson coeffi- 
cients of penguin operators in the effective Hamiltonian, q ~ 0.04 (i = 3, 4, 5, 6) (or 
a QCD loop factor), and a suppression by the Okubo-Zweig-Iizuka (OZI) rule [13]. 
The parameter £ is expected to be somewhat larger than the product of these three 
factors due to color-suppression of the dominant tree amplitude which normalizes £. 
Thus, with £ ~ 10~ 3 , it has been commonly accepted that the measurement of sin 2f3 
in B° — > J/ipK S ( L } may involve only a very small uncertainty at a level of 10 -3 , 
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or at most a fraction of a percent [6]. This estimate was supported by calculations 
of £ using QCD factorization [H] and perturbative QCD [T5]. These perturbative 
calculations are based on the absorptive part of the w-quark loop, assuming a suffi- 
ciently large momentum transfer in the loop, and applying rather crude methods for 
evaluating the four quark operator matrix element (J/ifjK°\(cr a c)v{bT a s)v-A\B ). 
The absorptive part of the quark loop was proposed thirty years ago as a mechanism 
producing a strong phase leading to CP violation in charged B decays [T6] . 

The introduction of a small penguin amplitude, carrying a weak phase 7 and a 
strong phase 5 relative to the dominant tree amplitude, affects the time-dependent 
CP asymmetry in B° — > J/ipKs, 

A M - H ~ -> r (A A^Q'/A * m 

Anpit) = n = — C cos Ami + b sm( Ami) , 1 

T(B° -> J/i/jKs) + T(B° - J/iPK s ) 

where, dropping terms quadratic in £ [6], 

C = — 2£ sin 5 sin 7 , AS = S — sin 2/? = 2£ cos 2/3 cos £ sin 7 . (2) 

In the limit £ = one has C = 0, A = 0. Current measurements of the two asymme- 
tries [31 HI [T7] (based on all charmonium decays), 

C(J/ijK°) = 0.005 ±0.019 , S( J/ipK°) = 0.671 ± 0.024 , (3) 

involve experimental errors at a level of ±0.02. This error is considerably larger than 
the theoretical uncertainty introduced by the above estimate of the parameter £. 

The estimate £ ~ 10~ 3 has been questioned and challenged in Refs. [H] and [19] . 
arguing that the w-quark penguin amplitude in B° — ► J/ipK° may be enhanced by 
long distance rescattering effects from intermediate 5 = 1 charmless states to J/ipK°. 
A sizable enhancement of a penguin amplitude beyond a perturbative calculation, ar- 
gued to be due to a large "charming penguin" contribution [20J, has been observed 
in B — > Kn. It was therefore argued [TSJ [TH] that similar nonperturbative rescat- 
tering effects of intermediate charmless states may enhance £, leading to a hadronic 
uncertainty in £ at a level of several percent. 

An uncertainty in £ at this level implies theoretical uncertainties in the asym- 
metries C(J/ipK°) and S(J/tpK°) which are comparable to or even larger than the 
current experimental errors in these measurements. It was pointed out in Refs. [T8lfT9] 
that CP asymmetries C and AS in B° — ► J/ipir , proportional to £ in the flavor SU(3) 
limit, are enhanced by a factor (1 - A 2 )/A 2 = 18.6 (A = 0.2257 [12]) relative to C and 
AS in B° — > J/ipKs- Thus, it was suggested to study £ in B° — > J/tpn . Unfortu- 
nately, the decay rate for this process is suppressed by 2A 2 /(1 — A 2 ) relative to that 
of B° — > J/tpKs- Consequently one expects errors in the B° — ► J/ipn asymmetries 
to be correspondingly larger than in the B° — ► J/ipKs asymmetries. Indeed, current 
measurements [13 M, [22] , C(J/^vr°) = -0.10 ± 0.13, S{J/^°) = -0.93 ±0.15, are 
not sufficiently accurate for providing useful information about £. Values of £ as large 
as a few percent cannot be ruled out by the two asymmetries. 

The purpose of this Letter is to calculate upper bounds on £ in B° — > J j ipK° from 
long distance rescattering effects mediated by charmless intermediate states. This 
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provides a re-evaluation of the contribution of the absorptive part associated with 
cutting the w-quark penguin loop [T4J [T5J [T6J by computing explicitly contributions of 
charmless intermediate states. Using the rich amount of data for numerous charmless 
B meson decays obtained in experiments at e + e~ B factories, we will show that the 
estimate £ ~ 10 -3 is much more reasonable than values of £ at a level of a few percent. 



II Upper bounds on rescattering in B° —> J/i/jK° 

We write the S matrix in terms of So, which includes strong and electromagnetic 
interactions, and T, taken to be Hermitian, which corresponds to the effective weak 
Hamiltonian at a low energy scale, 

S = S + iT . (4) 

Unitarity of the S matrix S^S = 1 implies to first order in T, 

T = S TS . (5) 

Taking matrix elements of the two sides between a B meson state and a final decay 
state /o, and inserting a complete set of intermediate states /, one has 

(fo\T\B) = V,{f \S \f){f\T\B), (6) 

where we used the fact that B is an eigenstate of Sq. The matrix elements (fo\T\B) 
and (f\T\B) are weak decay amplitudes, often denoted A(B — > f ) and A(B — > f), 
while (/o|So|/) represents a rescattering amplitude from / to Jo- 
Let us first consider the matrix element between B = B° and /o = J/ ipK° for the 
effective \AS\ = 1, AC = operator T u involving a CKM factor V* b V us , 

(J/^K°\T U \B°) = Z f (J/^K°\S \f)(f\T u \B ) . (7) 

Because B is a spinless particle the J /if) and K° are in a P-wave. Consequently the 
states / are all J = 0, P = — 1 S = 1 states. Since we are replacing the absorptive 
part associated with cutting the u-quark penguin loop by contributions of physical 
intermediate states, we consider only charmless states. This includes a long list of 
states, such as / = K*+n~, p~K + , K*°ir°, p°K°, uK°, K*°7], K*°rf, and K* Q + { 1430) tt~, 
but excludes K + tt~, K°tt° in an S-wave state and K* + p~,K*°p° in S and D waves 
which have P = +1. 

Parity and time-reversal symmetry of So imply a reciprocal detailed-balance rela- 
tion (we are assuming a single polarization state because J = 0), 

KJ/^ |So|/)| = |(/|So|J/^°)| . (8) 

Upper bounds on matrix elements \{f\So\J/ipK°)\ for each of the above states, / = 
K* + n~, _K"q + (1430)7t~, may be obtained using the following considerations. 
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We apply Eq. ([6]) to matrix elements between B° and the above charmless final 
states / for the effective |Ai>| = 1,AC = operator T c involving a CKM factor 

v; b v CSJ 

(f\T c \B°) = J: k (f\So\k)(k\T c \B ) . (9) 

The left-hand-side is dominated by a penguin amplitude which obtains a sizable 
"charming penguin" contribution [20J. Assuming that a single intermediate state 
k = D*~D+ can at most saturate the sum on the right-hand-side, one has 

\(f\S \D*-Dj)\\(D*-D+\T c \B°)\ < \(f\T c \B°)\ . (10) 

On the left-hand-side of (11 01) we can replace the D*~D+ state by J/i/jK°. One ex- 
pects \(f\So\J/tpK°)\ < \(f\So\D*~D+)\ because the first amplitude is OZI-suppressed. 
To calculate the ratio of B° decay amplitudes into D*~D+ and J/tpK° we use the 
expression for decay rates, 

where p*j is the momentum of one of the two outgoing particles in the B° rest frame. 
The measured branching ratios and the corresponding momenta are [12] 



B(B° -> D*~D+) = (8.3 ± 1.1) x 10" 3 , p* D + = 1735 MeV/c , 

B(B° -> J/ipK ) = (8.71 ± 0.32) x 10~ 4 , p* R0 = 1683 MeV/c . (12) 

This leads to \(D*-D+\T c \B )\/\(J/^jK \T c \B )\ = 3.04 ± 0.21. Using the central 
value, Eq. (flOj) may be replaced by 

| ( /|5 | W |<i ^)g>l 0)| . (13) 

We denote 

'- \(f\T'\B°)\ ] 1 ' 

and note that (f\T c \B°) is approximately the total B° decay amplitude into /, 
(f\T\B°); similarly (J/tpK°\T c \B ) w (J/^K°\T\B°). Combining Eqs. © and (fT5|), 
one then obtains the following upper bound on each of the terms contributing to the 
sum in (J7J), normalized by the B° decay amplitude into J/ipK°: 

_ \(J/^K Q \S \f)(f\T"\B»}\ I ( \(f\T\B°)\ V 
?/ \(J/4jK°\T\B°)\ 3 f \\(J/i>K°\T\B )\) ' 1 ' 

This upper bound is a central result in our analysis. It should be considered a 
strong inequality (in which the factor 1/3 may be replaced by 1/10) because it is 
based on a conservative inequality fflOl) and on presumably strong OZI-suppression 
of \{f\S \J/ipK°)\ relative to \(f\S \D*' Df)\. 
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Table I: Expressions for matrix elements (f\T u \B°) and (f\T c \B°) in B° — > VP in 
terms of graphical amplitudes. 



Final state 


{f\T u \B } 


(f\T c \B°) 




—T' 


-P'p 


p-K + 


-n 


—P' 




-C'y/^/2 


(P'p - Pew,v)/V2 


p°K° 


-C' P /V2 


(Pv - P'ew,p)/V2 


coK° 


C' P /V2 


(P{s + 2S P + lP^ W;P )/y/2 


K*°rj 


-cyVs 


{Py — P'p — S' v — \P'ew,v)I "^3 




c v /VE 


(2Py + P'p + AS'y - lP' EW y)/VG 



III Numerical upper bounds on rescattering 

We now study numerical bounds on rescattering parameters £f for numerous interme- 
diate states / in B° — > / — > J/tpK°. We start by discussing S = 1 charmless states 
/ = VP consisting of pairs of vector and pseudoscalar mesons. We use 

( \(f\T\B°)\ V_ B(B° — » f) (f KQ \ 
\\(J/^K \T\B°)\) B(B° -> J/ipK°) \p} ) ' 1 ' 

Values for the parameter ry, the ratio of two amplitudes in B —>■ VP involving 
CKM factors V* b V us and V* b V cs , are extracted from a study applying broken flavor 
SU(3) to these decays and decays into corresponding S = charmless states [23]. In 
the language of Ref. [21], matrix elements (f\T u \B°) involve combinations of graphical 
amplitudes representing color-favored and color-suppressed tree amplitudes T V , P \ and 
C'yrpy while (/|T C |_B°) involve penguin amplitudes Py/p), singlet penguin amplitudes 
S'y/ P \ (corresponding to SU(3) singlet mesons in the final state), and electroweak 
penguin amplitudes P'ew,v(p)- ^he subscript V or P denotes the final-state meson 
(vector or pseudoscalar) incorporating the spectator quark. We are neglecting color- 
suppressed electroweak penguin contributions. SU(3) breaking is included in T v and 
T'p in terms of ratios of pseudoscalar and vector meson decay constants, fx/ fit an d 
fx* I fp, respectively. Expressions for the above matrix elements are given in Table I 
for B° decays into seven VP states. Note that that while (f\T u \B°) involves color- 
allowed tree amplitudes T' P and T v for / = K* + ix~ and / = p~K + , it is governed by 
color-suppressed amplitudes C' P and C' v for all other final states. Consequently the 
values of rj in the first two processes are expected to be considerably larger than in 
the others. 

We calculate numerical values for rj using entries in the third column of Table V 
in Ref. [23], updating some values by fitting to more recent measurements of B(B° — > 
K*°tt°), B(B° -> K*°r]) and B(B° -> K*°r]'). Branching ratios for the seven B° -> VP 
decays [T7], corresponding center-of-mass momenta p* [T2], and values of ry are used 
to calculate from Eqs. (1X51) and (IT5|) upper bounds on £f for the these intermediate 
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Table II: Branching ratios [T7|, center-of-mass momenta [12] , parameters rj and upper 
bounds on £y for seven charmless intermediate VP states. 



Mode 


B 


P* 


r f 


Upper bound on 


/ 


(io- 6 ) 


(MeV) 




(IO" 4 ) 




10.3±1.1 


2563 


0.31±0.03 


7.9±1.1 


p-K+ 


8.6±1.0 


2559 


0.26±0.03 


5.6±1.0 


K*°ir° 


2.4±0.7 


2562 


0.09±0.04 


0.6±0.3 


p°K° 


5.4±1.0 


2558 


0.04±0.03 


0.5±0.4 


uK° 


5.0±0.6 


2557 


0.04±0.03 


0.5±0.4 


K*°T] 


15.9±1.0 


2534 


0.04±0.02 


1.6±0.7 


K*°7]' 


3.8±1.2 


2471 


0.08±0.04 


0.8±0.4 



VP states. Input values and resulting upper bounds on are summarized in Table 
HI The largest upper bounds, £/ < (7.9 ± 1.1) x 10~ 4 and £ f < (5.6 ± 1.0) x 10~ 4 , 
are obtained for / = K* + tt~ and / = p~K + , respectively. Much smaller values, at a 
level of 10 -4 , are calculated for all other VP states. 

Because rescattering effects of the form B° — > / — > J/ttK° increase with B(B° — > 
/), we search for charmless intermediate states / for which this branching ratio is 
particularly large. We note that the three-body decay mode B° f^°7r + 7r~, with 
B = (44.8±2.6) x 1(T 6 is dominated by the quasi-two-body decay B° -> K*(1430) + 7r~, 
involving a scalar and a pseudoscalar meson in a P = —1 S- wave |12j : 

B[B° -> K*(1430)+7r-] = (50^) x IO" 6 . (17) 

The fact that this branching ratio seems to exceed that for the three-body final 
state indicates strong destructive interference with other amplitudes including B° — > 
K* + ir~, p°K°, /o(980)ff°, and a non-resonant amplitude [25] . 

In order to evaluate an upper bound for £/ based on the w-quark amplitude's 
contribution to the Kq (1430) + 7t _ intermediate state, we must obtain an estimate 
of the value of r/ for this state. This quantity (the subscript P denotes the final 
pseudoscalar meson incorporating the spectator quark), 



\T' P (B° -> KZ+TT- 
\P' P {B Q -f K£ + ii- 



is the ratio of the it-quark tree amplitude and the c-quark penguin amplitude in 
B° — ► i^Q(1430) + 7r _ . While the latter amplitude dominates this process, the former 
may be estimated to a good approximation assuming factorization. A similar situation 
occurs in B° — > K + tt~, where the ratio of tree and penguin amplitudes has been 
determined within a global flavor SU(3) fit to all B — > Kit and B — » tttt decays 



\T'(B° -> K+7r-)\ 0.281(16.1 ±2.0) , . 

; = 0.094 ±0.012 . (19) 



\P'(B° -> K+w-)\ 48.2 ±1.0 
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The ratio of the two penguin amplitudes dominating B° — > K + ix and B c 
i^Q(1430) + 7r~ is obtained from the corresponding partial rates, 



\P'(B° -> K+n 



' " ; = 0.60 ±0.05 . (20) 



\P' P {B Q -> ^o + 7r-)| ^ -> i^o+Tr-) p*(X 

Here we have used with [12J -6(5° -> K+tt-) = (19.4±0.6)x 10~ 6 , p*(.K"+ #5+) = 
(2615, 2445) MeV/c. 

In the factorization approximation the tree amplitudes T'(B° —>■ K + tt~) and 
T'p(B° — > K* + n~) involve respectively the K and K$ decay constants, fx and fx*, 
and the B to 7r form factors at g 2 = m|- and m K * which are assumed to be approxi- 
mately equal. Thus, 

\T' P {B -f K*+tt-)\ ^f KS , , 

\T'(B° ^ K+tt-)\ ~ f K ■ 1 ) 

Note that the scalar Kq couples to the weak vector current through a coupling pro- 
portional to the Kq decay constant fx* which vanishes by G-parity in the SU(3) 
symmetry limit [271 [23 [29]. SU(3) breaking leads to a nonzero value, expected to 
be of order (m s — md) /Aqcb relative to the K meson decay constant. Theoretical 
calculations of fx* lead to values in the range 



fx* = 40 ± 6 MeV , (22) 



to be compared with fx = 155.5 ± 0.8 MeV 

Taking a product of the three factors in (fT9l . (120]) . and f[2Tj) . we find 

r f = 0.015 ± 0.003 , / = AT*(1430) + 7r . (23) 

An upper bound on for / = _K"q(1430) + 7t~ is then obtained using Eqs. (JT5|) and 
f fT6|) while taking into account correlated errors, 

. 1 B[B° -> K*(1430) + 7r-] p*(K°) 



3 1 B(B° -> J/^if ) p*(K *+) 
= (1.9 ± 0.4) x 10" 4 , / = fT *(i430) + 7r- . (24) 

Thus, in spite of the large branching ratio measured for this decay mode, this upper 
bound is about four times smaller than the largest value obtained for the correspond- 
ing VP state K* + TT~ in Table HH 

There are good prospects for replacing the theoretical estimate f[2"2l with an ex- 
perimentally determined value. The partial width for the decay r~ — > M~v, where 
M is a strange scalar or pseudoscalar meson, is 

r(r" - M-v) = Gl\V us \ 2 ^ { ^—^ . (25) 

167r m T 

With G F = 1.16637(1) x lO" 5 GeV" 2 , m T = 1.77684(17) GeV/c 2 , m K = 0.493677(16) 
GeV/c 2 , |K. I = 0.2255(19), and f K = 0.1555(8) GeV [12] (for the last two, see Ref. 
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[33!). this y ields a prediction r(r - -> tf _ z/) = (1.59 ± 0.03) x 10~ 14 GeV. The 
lifetime of the r is (290.6 ± 1.0) x 10~ 15 s |X2J , implying the prediction Bir~ — > 
K~v) = (7.02 ± 0.14) x 10 -3 , in satisfactory agreement with the experimental value 
P2] (6.95 ±0.23) x 10~ 3 . 

The prediction (125]) can be applied to the _K"q(1430) in the narrow-width approxi- 
mation, permitting one to obtain the ratio of scalar and pseudoscalar decay constants: 



Ik* _ m 2 T - m\ 
fx m? - m K * 

With the experimental upper limit [35] B\r -> K *(1430)z/] < 5 x 10~ 4 (95% c.l.) and 
using the central value of m[iT^(1430)] = 1425 ± 50 MeV/c 2 [12], this ratio is less 
than 0.69, entailing fx* < 107.5 MeV if the predicted value of B{r~ — > K~v) is used. 
The large width T[Kq (1430)] = 270 ± 80 MeV leads to a small positive correction of 
1.127 to the predicted partial width for r — > /<"o(1430)z/, reducing this upper bound 
slightly to fx* < 101 MeV. This is about a factor of 2.5 larger than the theoretical 
estimates summarized in Eq. ff22l . Using those estimates and including the finite- 
width correction, we predict B[t — > fQ(1430)z/] = (7.8 ± 2.3) x 10~ 5 . As stressed in 
Ref. [28], this should be accessible in present experiments. 

The intermediate state / = Kq (1430)°7r° is fed by a color- suppressed w-quark tree 
amplitude. Its penguin-dominated branching ratio is expected to be about half of 
that measured for B° — > Kq (1430) + 7r~. Consequently the upper bound on for 
/ = ifo(1430)°7r° is considerably smaller than (1241) . 

IV Conclusion 

We have calculated upper bounds on contributions to the doubly-CKM-suppressed 
parameter £ from rescattering B° — > / — > J/ipK° through charmless S = 1 inter- 
mediate states /. We have derived in Eq. (|T5l) a general conservative upper bound 
on £f which increases with B(B° — > /) and with the ratio r/ of tree and penguin 
amplitudes in B° — > f. The actual upper bound may involve a factor 1/10 instead of 
1/3 because OZI suppression in / — > J/ipK° has not been included in f|T5|) . 

The highest upper bounds on £f, somewhat below 10~ 3 , were obtained for / = 
K* + tt~ and p~K + , while other intermediate states with neutral vector and pseu- 
doscalar mesons involve much smaller rescattering contributions. This applies also to 
the state Kq (1430) + 7r _ , which has the largest quasi-two-body decay branching ratio 
measured so far in B° decays, B[B° -> fsT *(1430) + 7r-] = (50±|) x 10~ 6 . We noted de- 
structive interference between B° — > Kq (1430) + 7r~ and other modes contributing to 
B° — > K°7t + 7t~ . This indicates potential destructive interference between rescattering 
contributions to £ of these intermediate states. 

One may wonder whether larger contributions to £ may originate in charmless 
intermediate states with multiplicity larger than three. Only two 5 = 1 charmless 
branching ratios comparable to that of B° —>■ Kq (1430) + 7r~ have been measured [36], 
B(B° -> K*°tt + ti-) = (54.5 ±5.2) x 10" 6 , B(B° -> K*°K + K~) = (27.5 ±2.6) x 10" 6 . 



B(t - 
~B\F 



. Ku) 



(26) 
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These involve quasi-three-body decays leading to four particles in the final state. 
P = — 1 projections of these states, with smaller branching ratios, may rescatter into 
J/i/jK°. Although it is difficult to calculate or measure the tree-to-penguin ratio 
rf for these states, we expect it to be no more than 0.1. The isospin relation [37] 
T(B° -> K*°tt + tt-) = T(B + -> K*+tt + ii-) which holds within 1.5a [T7J is consistent 
with rf = in these processes. We do not anticipate constructive interference between 
rescattering contributions of the intermediate states K* 7i + 7T^ and K*°K + K~ and the 
above calculated contributions of K* + n~ and p~K + which are probably larger. 

Rescattering from intermediate states with two vector mesons in a P = —1 P- 
wave state, including K* + p~,K*°p° and K* <p, involve branching ratios considerably 
smaller than those of B° — > K*°ir + n~ and B° — »■ K*°K + K~ [12] and very small 
w-quark tree amplitudes. For instance, the tree amplitude in B° — > K* + p~ is related 
by flavor SU(3) to the amplitude dominating B° — ► p + p~ ■ Approximately 100% 
longitudinal polarization has been measured in this process [17J, corresponding to a 
combination of S and D waves but no P wave. This implies a negligible w-quark P 
wave amplitude in B° — > K* + p~. Similarly, the tree amplitude in B° —>■ K*°p° is 
color-suppressed, while the one in B° — > K* <p is both color and OZI-suppressed. The 
contributions of two vector meson intermediate states to £ are therefore negligible. 

Thus, we expect a value of £ which is at most a few times 10 -3 , in agreement with 
an early estimate [6] and in contrast to a suggestion for an order of magnitude larger 
enhancement of £ by long distance effects [TH [19] . 
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